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Abstract 



The dependence of the nuclear force on standard model parameters plays 
an important role in bounding time and space variations of fundamental cou- 
plings over cosmological time scales. We discuss the quark-mass dependence of 
deuteron and di-neutron binding in a systematic chiral expansion. The lead- 
ing quark-mass dependence of the nuclear force arises from one-pion exchange 
and from local quark-mass dependent four-nucleon operators with coefficients 
that are presently unknown. By varying these coefficients while leaving nu- 
clear observables at the physical values of the quark masses invariant, we find 
scenarios where two-nucleon physics depends both weakly and strongly on 
the quark masses. While the determination of these coefficients is an exciting 
future opportunity for lattice QCD, we conclude that, at present, bounds on 
time and space variations of fundamental parameters from the two-nucleon 
sector are much weaker than previously claimed. This brings into question 
the reliability of coupling-constant bounds derived from more complex nuclei 
and nuclear processes. 



I. INTRODUCTION 



The recent observation suggesting that the fine-structure constant was smaller in the 
past in than it is today has led to renewed interest in the idea of using time (and space) 
variation of fundamental parameters as a probe of high-energy physics. Based on the prin- 
ciple that all that is not forbidden is mandatory, time variation in the electromagnetic force 
is perhaps not surprising. Moreover, one would expect that the other forces of nature vary 
with time as well. This idea is quantifiable by assuming a suitable grand unified theory 
(GUT) at a scale Mqut or a particular brane-world scenario The assumption that the 

short- distance couplings of the theory are related at all times leads to relations between the 
variation of the strength of the electromagnetic interaction, aem, and the variation of other 
standard model parameters, such as the light quark masses, m^. As a change in and 
aem will naively lead to a change in the positions of nuclear energy levels, special interest 
has been paid to the light-element abundances predicted by big bang nucleosynthesis (BBN) 
and also to the abundance of isotopes produced by the Oklo "natural reactor" in the hope 
that these abundances can be used to constrain high-energy physics [p|,^,|5|- p^] . 

In this work we will critically analyze the two-nucleon sector using an effective field theory 
(EFT) that respects the approximate SU{2)l ® SU{2)r chiral symmetry of QCD and has 
consistent power-counting This theory allows for a systematic study, consistent with 

QCD, of the bound-state in the ^Si — ^Di coupled-channels, the deuteron, and the scattering 
amplitude in the ^5*0 channel. This is the first study of the mg-dependence of these quantities 
that is complete at next-to-leading order (NLO) in a consistent EFT. Unfortunately, due to 
the current lack of understanding of the low-energy behavior of QCD, only weak bounds, if 
any, can be placed on the time-dependence of fundamental couplings from the two-nucleon 
sector. We will explicitly demonstrate this by finding plausible values for couplings in the 
EFT for which the deuteron binding energy varies relatively little, and for which the di- 
neutron system in the ^5*0 channel is unbound, over a relatively wide range of quark masses. 

The deuteron is an interesting object to study for a variety of reasons. From a fundamen- 
tal point of view it is quite intriguing as it is bound by only = 2.224644 ±0.000034 MeV, 
which is much smaller than the typical scale of strong interactions, Aqcd- Clearly, in order 
to arrive at such a small binding energy, fine-tunings are involved f\ and consequently naive 
dimensional analysis (NDA) as applied to the single-nucleon sector or the mesonic sector, 
which we will call naive naive dimensional analysis (N^DA), on this system is doomed to 
fail. Unfortunately, all previous attempts to extract bounds on the variations of fundamen- 
tal couplings from the deuteron, and in general the two-nucleon sector, have implemented 
N^DA. From a more phenomenological standpoint, the smallness of plays a key role in 
the synthesis of light elements in BBN. The impressive agreement between the predictions 
of BBN and observation suggests that new physics that would have significantly modi- 
fied the deuteron at the time of BBN is absent. The ^5*0 channel is quite similar to the 
^5*1 — ^Di coupled-channels in one important way, its scattering length is unnaturally large, 



*In the pionless EFT one finds that both the and the ^Si — ^Z?i coupled-channels are very 
close to an unstable IR fixed point of the renormalization group. At this fixed point, the scattering 
lengths are infinite [26,2^. 
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a^^So) — —23.714 ± 0.013 fm. While there is no bound state in this channel for the physical 
values of the light-quark masses, a small increase in the strength of the nuclear force would 
bind two nucleons in this channel. The existence of a bound state in this channel, e.g. a di- 
neutron, nn, in the nucleosynthesis epoch would be quite profound and could substantially 
modify the predictions of BBN. 



II. EFT FOR TWO NUCLEONS 

During the last decade there has been a significant effort to construct an EFT to describe 
nuclear physics. While it is straightforward to write down all possible terms in the effective 
Lagrange density for two or more nucleons, arriving at the correct power-counting proved to 
be a difficult task. Weinberg's (W) original proposal [0] for an EFT describing multi-nucleon 



systems was to determine the nucleon-nucleon (NN) potentials using the organizational 
principles of the well-established EFT's describing the meson-sector and single-nucleon sector 
(chiral perturbation theory, xPT), and then to insert these potentials into the Schrodinger 
equation to solve for NN wavefunctions. Observables are computed as matrix elements 
of operators between these wavefunctions. W power-counting has been extensively and 
successfully developed during the past decade to study processes in the few-nucleon systems. 
This method is intrinsically numerical and is similar in spirit to traditional nuclear-physics 
potential theory. Unfortunately, there are formal inconsistencies in W power-counting p5| , 



in particular, divergences that arise at leading order (LO) in the chiral expansion cannot 
be absorbed by the LO operators. Problems persists at all orders in the chiral expansion, 
and the correspondence between divergences and counterterms appears to be lost, leading 
to uncontrolled errors in the predictions for observables. This formal issue was partially 
resolved by Kaplan, Savage and Wise (KSW) who introduced a power-counting in which 
pions are treated perturbatively [| [^]. The NN phase-shifts and mixing angle in the ^5*0 and 



^5*1 — coupled-channels have been computed to next-to-next-to-leading order (N^LO) 
in the KSW expansion by Fleming, Mehen and Stewart (FMS) from which it can be 
concluded that the KSW expansion converges slowly in the ^5*0 channel and does not converge 
in the ^Si — ^Di coupled-channels. Therefore, neither W or KSW power-counting provide a 
complete description of nuclear interactions (for recent reviews see Ref p0|-[3^ 



The problems with W and KSW power-counting appear to have been resolved in the 



work by Bedaque, van Kolck and the authors which from this point on we will refer 
to as BBSvK. It was realized in FMS that the contributions to the amplitude that lead to 
non-convergence in the ^Si — ^Di coupled-channels persist in the chiral limit (it is the chiral 
limit of iterated one-pion-exchange (OPE) that is troublesome). Therefore, in BBSvK power- 
counting the scattering amplitude is an expansion about the chiral limit. This recovers KSW 
power-counting in the ^5*0 channel, where FMS found it to be slowly converging. However, 
in the ^5*1 — ^Di coupled-channels, the chiral limit has contributions from both local four- 



^ In KSW power-counting the momentum-independent four-nucleon operators are LO and are 
resummed to yield the LO scattering amplitude, while pions are NLO and treated in perturbation 
theory. 
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nucleoli operators and from the chiral limit of OPE. It is these two contributions that must 
be resummed using the Schrodinger equation to provide the LO scattering amplitude in the 
'^Si — '^Di coupled-channels |. 

The LO Lagrange density describing the single-nucleon sector and the pseudo-Goldstone 
bosons in two-flavor QCD is 



£ = -(/f ))2 Tr [ 1 + ATr [ m^S^ + m^S 

g L J L 

+ Art + W^/2Mf ) N + gfN^a ■ AN , 



(1) 



where /^°^ is the pion decay constant in the chiral limit, Mj^ is the nucleon mass in the 
chiral limit, g^^^ is the axial coupling constant in the chiral limit, E is the exponential of 
the pion fleld, and A is the axial-vector field. The properties of the nucleons and mesons 
have been studied in the first few orders of the chiral expansion (a good review can be found 
in Ref. [^). If one is interested in the mg-dependence of the nuclear force, as we are, one 
needs to have the chiral expansion for the nucleon mass, for the axial coupling and for the 
pion decay constant up to NLO. Each of these observables has been studied extensively, the 
results of which can be found in Refs. 
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where I4 = 4.4 ± 0.2 |p7| , ^ , ci ~ — 1 GeV~^ are mg-independent constants, and = 
139 MeV has been used to determine the constants in the chiral limit. We use qa = 1.25, 
Mm = {Mn+Mp)/2 and /,r = 135 MeV. We have retained only the leading chiral-logarithmic 
contribution to qa, and have chosen a renormalization scale of A = 500 MeV. This point 
requires discussion. Extraction of the counterterm relevant to the mg-dependence of qa at 
one-loop order presently yields A ~ 100 MeV (see Ref. and references therein). As this 
anomalously small value seemingly indicates a breakdown of the chiral expansion, we assume 
that this is a problem with the extraction and use a natural value of A = 500 MeV for this 
analysis. This, of course, introduces an uncertainty at NLO in the EFT calculation. 



A technical issue that remains to be solved is that while dimensional regularization can be used 
to regulate the channel, so far no one has managed to dimensionally regulate the "^Si — ^Di 
coupled-channels due to the presence of the component of the nucleon-nucleon potential. 

Therefore, in this work we will dimensionally regulate the ^Sq channel and use position-space 
square-well regularization in the ^Si — ^Di coupled-channels. It has been shown that a singular 
l/r" potential (n > 2), can be regulated and renormalized by a single momentum-independent 
square- well p5[ |. 
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The interactions between two nucleons arise from pion exchange, resulting from the 
Lagrange density in eq. (|l|), and also from local four-nucleon interactions, which for s-wave 
interactions result from a Lagrange density of the form. 



Ids2 (n^n) (iVtiv) Tr [M,+] - \dt2 (iVVw) [N^a'N) Tr [M,+] + ... , (3) 



2 " V / 2 
1 

2 
1 

~2 

where M.q+ = | ( i'^^qC + ^f^qi ); and = diag(m„,md). The ellipses denote operators 
involving two derivatives [|^, which are the same order in the power-counting as a single 
insertion of m^, and also higher- dimension operators. The rrig-independent coefficients Cj, Di 
and so forth, are to be determined from experimental data. An important point to notice is 
that one cannot separate the contributions from the Ci and the Di using NN scattering data 
alone, as both operators are momentum independent. However, the Di operators contribute 
to interactions between two nucleons and two or more pions, while the Ci operators only 
contribute to interactions between two nucleons. 



III. THE ^Sq channel AND THE DI-NEUTRON 

One might imagine that grave violence would be done to the predictions of BBN if a 
di-neutron or di-proton were to be stable during the nucleosynthesis epoch. Thus it is 
interesting to know just how much the couplings in the Lagrange density describing the 
interactions between two nucleons in the ^5*0 channel can change before a di-neutron state 
becomes stable 0. However, from the point of view of setting bounds on the variation of fun- 
damental couplings, what is important to know is if there is a set of "reasonable" couplings 
for which the scattering length remains unnaturally large and the di-neutron is unbound 
over a large range of pion masses. In this situation the pion mass can vary substantially, yet 
still not significantly modify the physics of two nucleons in the ^5*0 channel. Indeed, such a 
parameter set exists. 




FIG. 1. LO contribution to the scattering amplitude in the Sq channel. The small solid circle 
denotes an insertion of Cq. 



§In this initial study, we do not consider the di-proton system due to the coulomb potential that 
further destabilizes it. However, for a meaningful study of constraints from BBN in this channel, 
the possibility of nn, np and pp bound states would have to be considered. 
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FIG. 2. NLO contributions to the scattering amplitude in the ^Sq channel. The small solid 
circles denote an insertion of Cq or (7^. Dashed lines are pions and the large solid circle (square) 
corresponds to an insertion of C2 (-D2). 



The scattering amplitude for two nucleons in the ^5*0 channel has been computed out 
to N^LO p6| , |29[| , but for our purposes it will be sufficient to work with the amplitude at 



NLO |26|. The scattering amplitude thus has an expansion of the form 

A = A-1+ Ao + Ai + . . . 



(4) 



where An is of order Q" in the dual and momentum expansions. The LO amplitude 
arises from the diagrams in Fig. ffl and is given by 



^-1 



while the NLO amplitude arises from the diagrams in Fig. ^ and is given by the sum of 
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where we have chosen to work in the isospin hmit, = ma, and we have turned off the elec- 
tromagnetic interaction. These amplitudes are manifestly renormalization-scale independent 
order-by-order in the EFT expansion. The coefficients Cq, D2 and C2 are the combinations 
of couplings from the Lagrange density in eq. (|^) appropriate for the ^5*0 channel. Cq corre- 
sponds to a four-nucleon operator that is independent of rriq and momentum, D2 corresponds 
to a four-nucleon operator that has a single insertion of and no derivatives, and C2 cor- 
responds to a four-nucleon operator that is independent of rrig and has two derivatives. The 
quantity fi is the renormalization scale, and we have used dimensional regularization and 
the power-divergence subtraction procedure (PDS) |^ to renormalize the theory; p is the 
magnitude of the three momentum of each nucleon in the center-of-mass frame. In addition, 
we have used the LO relation between the quark masses and the square of the pion mass 
so that the amplitude at NLO is written entirely in terms of m^r. It is important to note 
that the operator with coefficient D2 is required at NLO. It is this operator that absorbs 
the scale dependence of ^o'^^ 

From the NLO amplitude it is easy to construct pcot6''^^°\ which has a well-behaved 
power-series expansion for p < and thus we can determine a linear combination of 

Cq and D2 in terms of the scattering length and C2 in terms of the effective range. 



2.73 ± 0.03 fm, at the physical value of the pion mass. Once these parameters are 



fixed, the scattering length, effective range and phase-shift can be determined as a function 
of the pion mass. The scattering length is given by 
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Given that it is only a combination of Cq and D2 that can be 
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FIG. 3. The scattering length in the ^Sq channel (in fm's) as a function of the pion mass, 
for the couplings given in eq. (P). As the scattering length is negative over this entire region, the 
di-neutron is unbound. 



fixed from NN scattering, and at present D2 cannot be isolated with pionic processes, unique 
values for D2 or Cq that are determined experimentally do not exist. Thus, one must explore 
a "reasonable" range of values for D2 and Cq in order to gain some understanding of how 
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the ^5*0 channel is modified away from the physical values of the quark masses. It is quite 
straightforward to identify a set of couplings, 

Co(m™YS) = -4.09 fm^ , 1^2 (m™^^) = 0.50 fm^ , Cslm™^^) = 4.94 fm^ , (8) 

for which the scattering length in the ^5*0 channel remains unnaturally large and negative 
over a wide range of pion masses, as shown in Fig. ^, and hence a di-neutron is not stable in 
this region We have extracted the Cj and D2 using eq. (0) and an analogous equation for 
the effective range, with qa-, Mn and /vr set to their chiral limit values in eq. (^. In order to 
determine the impact of a varying pion mass on nuclei more complex than the deuteron, it is 
useful to know the behavior of the phase-shift as a function of the pion mass over a relatively 
wide range of momentum. In Fig. ^ we show the ^Sq phase-shift for m,r = 60 MeV, m^^^^ 
and 180 MeV, for the couphngs in eq. (H). One notices that the NLO phase-shift does not 




FIG. 4. The phase-shift, (J^^"^"-*, as a function of momentum, |p|, for pion masses of 
= 60 MeV (dashed), ml^^^ (dotted) and 180 MeV (dot-dashed), for the couplings in eq. (§). 



The solid curve corresponds to the results of the Nijmegen partial-wave analysis |41|. 



reproduce the results of the Nijmegen partial- wave analysis |iT| above ~ 50 MeV, even at 



the physical value of the pion mass |^. This is because we have fixed the parameters, Cq, 
D2 and C2 to the scattering length and effective range. However, this means that the dual 
chiral and momentum expansion of the EFT at NLO has been matched to the momentum 
expansion of pcot^^'^^^. Consequently, both the scattering length and the effective range 
have their own chiral expansions that we are truncating at NLO. Thus we do not expect to 
reproduce the phase-shift exactly, but should be perturbatively close, as can be seen to be 
the case in Fig. |^. Rather than compare the phase-shifts with the Nijmegen partial-wave 
analysis El|, it is perhaps more informative to compare the phase-shifts for = 60 MeV 



**For convenience, we plot observables as a function of 771,7^ rather than oc ruq, which would 
be more natural from the standpoint of QCD. 



^For a detailed discussion of the ^5*0 channel description with EFT, see Refs. [ p6| , pS| , 42 1 . 
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and 180 MeV with those at m,r = . Thus we have identified a set of NLO couphngs 

for which the ^5*0 channel is quite insensitive to moderate changes in the pion mass. 

It is interesting to ask what one might expect if nature has chosen values for the couplings 
other than those in eq. (|). Choosing some arbitrary values for the couphngs, but ones that 
still respect NDA in the two-nucleon sector [p^ , 

Co(m™YS) = -4.00 fm^ , Dal^n™^^) = 0.31 fm^ , Calm™^^) = 4.57 fm^ , (9) 
we show the scattering length as a function of the pion mass in Fig. |^. One can see that there 




FIG. 5. Properties of the ^Sq channel for the couplings in eq. @. The left panel shows the 
scattering length as a function of the pion mass. For these couplings, there is a window where 
the di-neutron is stable. The right panel shows the phase-shifts for ttt-tt = 60 MeV, m™^^ and 
180 MeV. 



are couplings for which the scattering length becomes positive, indicating the presence of a 
di-neutron that is stable with respect to the strong interaction, for relatively small variations 
in the pion mass. Perhaps the unjustified choice of NN interaction arising from N^DA used 
in previous works on this subject (e.g. Ref |T7| , p!8| , pT| ) is in some way comparable to such an 
ad-hoc choice of couplings in the EFT. 

It is clear that sufficiently little is known about the m^-dependence of the ^5*0 channel 
that placing bounds on the time-variation of from nuclear processes sensitive to this 
channel is not possible at this point in time. 



IV. THE DEUTERON BINDING ENERGY 

As we have already discussed, the deuteron plays a central role in BBN and the produc- 
tion of light elements. If the deuteron was deeply bound or unbound in the nucleosynthesis 
epoch then the abundances of the light elements would look radically different from the 
predictions of BBN and from those observed in nature. Therefore, it is possible that limits 
can be placed upon the time-variation of fundamental couplings if the dependence of i?d, 
and relevant nuclear cross-sections on these couplings is known. First realistic attempts 



to understand Bf^{mq) as a function of were made in Refs P7|j43[| . One motivation for 
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BBSvK was to understand just what was possible for the mg-dependence of Bd{mq), and in 
particular to determine if the deuteron was bound or unbound in the chiral limit. This issue 
could not be resolved due to the lack of information about the D2 operators that contribute 
to the ^Si — ^Di coupled-channels. A second motivation was to understand what was required 



to make any sort of contact with quenched lattice QCD calculations of the scattering 
lengths in the ^5*1 channel. 






FIG. 6. Lippmann-Schwinger equation for the LO contribution to the scattering amplitude 
(large solid rectangle) in the ^Si — ^Di coupled-channels. The small solid circles denote an insertion 
of Co or QA- The "o" appearing below the OPE diagram implies keeping only the chiral limit. 



In the '^5'i — ^Di coupled-channels, OPE generates both central and tensor potentials, 

r(7r)/ _ N _ 2 ^ T/('^)/' . ^ _ (1 "^''^1" 



Vc ij-m^) = -aT,m^ — - — , {r]m-„) = -a^ — - — ( ^ + — + ) , (10) 



where a^, = g\{l — 2ml.dis/gAY/ (Stt/^). The constant dis is uncertain |^ but as our results 
are quite insensitive to this quantity, we set it to zero in what follows. At LO in BBSvK 
counting we should only keep the chiral limit of these potentials and the Co operator |^ that 
contributes to the ^Si — ^Di coupled-channels, as shown in Fig. ^ Deviations from the chiral 
limit are inserted in perturbation theory. We choose to keep the full potentials in our NLO 
calculation, meaning that we also include part of the N^LO calculation, which we find does 
not modify the divergence structure. In addition, at NLO in BBSvK counting there is a 
contribution from the chiral limit of two-pion exchange (TPE) and from an insertion of C2 
and D21 as shown in Fig. |^. The TPE potential in coordinate space has been computed in 



Ref. [45|, and in the chiral limit is given by 



^...,(,0)^ 3(22.1 -J0.1-l) J^ . Vt\r,0) ^ - . (U) 

In order to regulate the potentials at the origin, we use a spatial square-well of radius 



R 46,27 , where the potential outside the square well is 



Min contrast with the ^So channel, in the ^Si — ^Di coupled-channels there is a local four-nucleon 
interaction from OPE that we have included in our definition of Cq. 
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2C-(^-^)-X-X 

o o 

FIG. 7. Chiral limit of the crossed TPE diagram, deviations from the chiral hmit of OPE, and 
the C2 (large solid circle) and D2 (large solid square) operators, all of which contribute at NLO 
in the — ^Di coupled-channels. The "o" appearing below a diagram implies keeping only the 
chiral limit. 



VL\T,m^) y_2^M^VTir;m^) -M^iVcir;m^)-2VTir;m^))-Q/r' J ' ^'^^ 
where 

Vc{r;m^) = vi''\r;m^) + V^''''\r;0) , VT{r;m^) = T^^^r; m^) + l^^^^r; 0) . (13) 
The potential inside the square well is 

V.fr-m fc^)- ('-^^(^'^o + ^dX + ^VcJ A 

^s{r,m^,k}-^ Q _M^^Vc, + Vn,ml + k'VaJ-6y) 

where Vcq, and are constant potentials corresponding to the renormalized local 
operators with coefficients Co, D2 and C2 in the ^Si — ^Di coupled-channels, respectively, 
and again we have used the LO relation between the pion mass and m^. It is important 
to recall that there is implicit mg-dependence in this potential due to qa, Mjy and /t,. For 
the chiral limit of the OPE potential, the NLO mg-dependence of these constants must be 
retained. However, for the deviations of the OPE potential from the chiral limit, and also 
for the chiral limit of the TPE potential it is sufficient (and in fact necessary to be able to 
renormalize the theory) to keep only the chiral limit of these constants. Defining to be 

where u{r) is the S-wave wavefunction and w{r) is the D-wave wavefunction, the regulated 
Schrodinger equation is 

^"(r) + + VL{r;m^) 9{r - R) + Vs{r;m^, k'^) e{R - r) ]<iJ{r) = 0. (16) 

As discussed in BBSvK, we can make an identification between the coefficients of local 
operators, Cj and Di, and the constant potentials of the square-wells that enter into eq. (0), 
Vc,; and Vd„ 

A 5''\r) ^^^H^ - V,^ e{R - r) . (17) 
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FIG. 8. The deuteron binding energy at NLO in the EFT as a function of the pion mass for 



the couphngs given in eq. (18). For this set of parameters, the deuteron is loosely bound over a 
wide range of pion masses. 



At the physical value of the pion mass, C2 and a linear combination of Co and D2 are 
fixed to reproduce the deuteron binding energy and the effective range in the ^Si-channel. 
Given that it is only a combination of Co and D2 that we determine, we have freedom to 
vary D2 and compensate this by a change in Cq, as in the ^Sq channel. We find couplings 
Co, D2 and C2, renormalized at the scale set by the radius of the square-well, R* = 0.45 fm, 

Co(i?*) = -6.17 fm^ , D2{R*) = 0.67 fm^ , C2{R*) = 0.75 fm\ (18) 

that are consistent with NDA and all the low-energy phase-shift data in the ^5*1 — ^Di 
coupled-channels, and for which the deuteron is loosely bound over a large range of pion 
masses, as shown in Fig. ^ 

The phase-shifts and mixing parameter for the ^Si — ^Di coupled channels are shown in 
Fig. |. As expected, the ^Si phase-shift falls more steeply as the deuteron becomes more 
loosely bound. Further, as the pion becomes lighter, the higher partial waves and mixing 
parameter are expected to increase due to the longer range of the potential and their relative 
insensitivity to short-distance physics. 

In BBSvK, the renormalization scale dependence of the theory was investigated ||27|| . 
The dependence on the width of the square-well, R, was found to be small, consistent with 
the EFT expectations. We find the same to be true at NLO. 

For purposes of contrast, we present the deuteron binding energy, phase-shifts and mixing 
parameter for an arbitrary set of couplings that respect NDA, 

Co(i?*) = -5.05 fm=^ , D2{R*) = -1.60 fm^ , C2{R*) = 0.75 fm^ . (19) 



^Ht is interesting to note that the chiral-limit value of the deuteron binding energy in the EFT 
coincides with that found in an analogous calculation with the Argonne V18 potential, when all 
coupling constants are frozen to their physical values, and only the explicit ttItt dependence in 
eq. dl^) is considered 
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FIG. 9. The phase-shifts, 6q and 62 for the ^Si channel and the ^Di channel and the mixing 
parameter Si as a function of momentum, |p|, for pion masses of m-n- = 60 MeV (dashed), m^^^^ 
(dotted) and 180 MeV (dot-dashed), for the couplings in eq. (jl^). The solid curve corresponds to 
the results of the Nijmegen partial-wave analysis @. Note that for 60 at ml^'^^, the NLO EFT 
calculation coincides with the partial-wave analysis to relatively high momenta. 
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FIG. 10. The deuteron binding energy at NLO in the EFT as a function of the pion mass for 



the couplings given in eq. (1£) 



For these choices, the variation of B^{mq) with respect to is quite rapid, as shown in 
Fig. |T^, and for m^r^ 90 MeV, the deuteron is unbound. The phase-shifts and mixing 
parameter resulting from the couplings in eq. ([T9|) are shown in Fig. |lT|, and as expected, 
they have a relatively strong dependence upon irtq. 

It is important to point out that while the deuteron binding energy curves shown in 
Fig. I and Fig. |10| yield unbound deuterons in the chiral limit, there exist parameter sets 
consistent with NDA for which the deuteron remains bound in the chiral limit. 



V. CONCLUSION 



We have explored the light-quark mass dependence of low-energy nucleon-nucleon inter- 
actions. The motivation for this work was to determine if, in fact, bounds could be set on 
the time- variation of fundamental couplings from nuclear processes, such as those occurring 
during big bang nucleosynthesis. We have demonstrated the existence of sets of strong in- 
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FIG. 11. The phase-shifts, 5o and 62 for the ^Si channel and the ^Di channel and the mixing 
parameter ei as a function of momentum, |p|, for pion masses of = 60 MeV (dashed), 139 MeV 
(dotted) and 180 MeV (dot-dashed), for the couplings in eq. ([l^). The solid curve corresponds to 



the results of the Nijmegen partial- wave analysis |41| 



teraction couplings in the low-energy effective field theory describing the nucleon-nucleon 
interaction that are consistent with all available data and with naive dimensional analysis for 
which the di-neutron remains unbound, and the deuteron remains loosely bound over a wide 
range of light quark masses. We do not mean to imply that these are the sets of couplings 
that nature has chosen, but rather that this scenario is not excluded at present. Thus, we 
conclude that bounds that have previously been set on the time-variation of fundamental 
couplings from processes in the two-nucleon sector are not rigorous and should be discarded. 

Our calculation does suffer from some limitations. First, we have not included electro- 
magnetism, and thus could not address the pp-system in the ^S'o-channel. However, we do 
not believe that its inclusion will modify the results we have presented in any significant 
way. The same we believe to be true of isospin breaking. Second, and perhaps the most 
important limitation, is that we have not included the strange quark in our discussions, 
and have worked with SU{2)l ® SU{2)r chiral symmetry. Variations in the strange (and 
charm, bottom and top) quark mass will manifest themselves as changes in the values of the 
coefficients in the Lagrange density, Ci and Di. At this point in time the tools do not exist 
to analyze this scenario. For instance, one might consider developing an EFT with three 
light-quark flavors. However, since nuclei are such finely-tuned systems, and SU (3)-breaking 
is not so small, it is likely that a perturbative treatment of nuclei using the approximate 
5'f/(3) L ® SU{3)r chiral symmetry of QCD would converge very slowly, if at all. 

From a conventional nuclear physics point of view, the mg-dependence of the nuclear 
potentials not only requires knowledge of the mg-dependence of one-pion exchange, but also 
of the rrig-dependence of the p, u, 0,... masses and their couplings to nucleons. This would 
appear to be an intractable problem. It is possible that some of the traditional treatments of 
the nucleon-nucleon interaction |5l,48-50|, and treatments of light nuclei could mimic these 
dependences by variations in their ad-hoc short-distance components of the nucleon-nucleon 
and higher-body potentials. The recently developed effective field theory tools organize this 
problem in a very simple way, and thus at NLO in the EFT expansion, only the Di are 
needed. In experimental processes examined to date that, in principle, allow one to separate 
the Di from the Ci, it is found that there are other amplitudes contributing to the process 
that dominate over the Di contributions, rendering such a separation exceedingly difficult. In 
the absence of an experimental determination of the Di, it would appear that the only viable 



13 



means for determination of this mg-dependence is through lattice QCD simulations |^ We 
find this to be very strong motivation to pursue a lattice QCD program focused on the 
two-nucleon sector. 
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